function [g] = g_function(p,data,wr)
% the moment function 
%   inputs: p - parameters
%           data - matrix of observed data
%           wr - estimate of reservation wage
%   outputs: g - matrix (m X N) where m is the number of moments
%                    E[g] = 0

% assign parameters
lambda = p(1);
delta = 

% assign data
u = data(:,1);
tu = 

% wage distributuion
F = 
f = 

% create moments
urate = 
uduration = 
eduration = 
ExpWage =
ExpWage2 = 

% create g vector
g = 

end

